Landau-Zener transitions in a multilevel system. An exact result. 
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We study the S-matrix for the transitions at an avoided crossing of several energy levels, which 
is a multilevel generalization of the Landau-Zener problem. We demonstrate that, by extending the 
Schrodinger evolution to complex time, one can obtain an exact answer for some of the transition 
amplitudes. Similar to the Landau-Zener case, our result covers both the adiabatic regime (slow 
evolution.) and the diabatic regime (fast evolution). The form of the exact transition amplitude 
coincides with that obtained in a sequential pairwise level crossing approximation, in accord with 
the conjecture of Brundobler and Elser [10] . 



Landau-Zener (LZ) transitions [1,2] between two en- 
ergy states at an avoided level crossing is one of the 
few exactly solvable problems of time-dependent quan- 
tum evolution. The LZ theory has found many applica- 
tions, and is used to describe a variety of physical sys- 
tems, ranging from atomic and molecular physics [3], to 
mesoscopic physics [4]. 

However, since most of the problems of interest involve 
more than two energy levels, with transitions between 
several levels happening simultaneously, it is desirable to 
extend the LZ theory to multilevel problems. A num- 
ber of interesting generalizations of the LZ problem have 
been proposed, in which several energy levels cross at the 
same time. In some cases, including the so-called bow-tie 
model [5-7], the high spin model [8], one can construct 
an analytic solution. The solutions [5-8] rely on a spe- 
cial form of the coupling Hamiltonian or on its symmetry, 
which preserves the integrability of the problem. One of 
the most efficient methods proposed to treat the general- 
ized LZ problems is the contour integration approach [9] . 

Much less is known about the general multilevel LZ 
problem, defined by a time-dependent Hamiltonian of the 
form 



H(t) = A + vt 



(1) 



were A, v are hermitian matrices nxnofa general form. 
In this problem, one is interested in the evolution 



(2) 



from large negative to large positive times. Since vt is 
the leading term in the Hamiltonian at large t, it is nat- 
ural to analyze the evolution in terms of an S-matrix in 
the basis of states that diagonalize v. 

One can construct an approximate solution in the 
'weak coupling' limit, when the off-diagonal elements of 
A are weak compared to the diagonal ones, Ay| <c |Ajj|, 
(i j). Since in this case the anticrossing gaps of the 
adiabatic levels (which sometimes are also called 'frozen 
levels') are much smaller than the level separation, the 
weak coupling regime can be analyzed in a sequential 
two-level LZ approximation [10], similar to the energy 
diffusion problem in mesoscopic systems [4]. 



Brundobler and Elser have noticed [10] that the for- 
mula for certain diagonal elements of the S-matrix ob- 
tained in such an approximation remains numerically 
very accurate even at strong coupling, when this approx- 
imation supposedly breaks down. This surprising obser- 
vation, based on a computer simulation of the problem 
(2),(1), has led Brundobler and Elser to a conjecture that 
the sequential LZ approximation can in some cases give 
an exact result. 

Below we show that this conjecture is indeed correct. 
We shall use the method of continuing the Schrodinger 
evolution (2) to complex time, and matching the t — > ±oo 
asymptotics of the evolving state ip(t). Our method gives 
exact results only for some of the transition amplitudes, 
corresponding to the transitions |l)t=_oo ~ * |l)t=+oo> 

™}t=-oo — > |»)t=+oo- 
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FIG. 1. Time evolution of the adiabatic (frozen) energy 
levels of the Hamiltonian (1). The transitions analyzed in 
this work are |l) t= -oo — » |l)t=+«> and \n)t=-oa — » |n) t=+00 . 

It is convenient to analyze the problem in the basis of 
states that diagonalize v: 



V-2 



\ 



(3) 



This basis is not the basis of adiabatic states, however, for 
t — * ±oo those states coincide with the adiabatic states, 



1 



up to permutation. It is also convenient to order the 
eigenvalues Vi so that 



(4) 



For simplicity, we assume that the eigenvalues are non- 
degenerate: v-i ^ Vj for i ^ j. 

To produce a full solution of the Landau-Zener prob- 
lem completely, one would have to find the transition am- 
plitudes between all different states, or, in other words, 
to find the scattering matrix S, relating the initial state 
of the system at t = — oo and the final state at t = +00: 



ip(t = +00) = Sip(t = -00) . 



(5) 



The Landau-Zener theory provides an analytic solution 
of this problem for n = 2. Presently, it is not known 
whether a generalization of this solution exists for n > 2. 

Here, I will pursue a less ambitious goal. I will show 
how one can obtain the two transition amplitudes Su 
and S nn . In our 'i/-ordered' basis notation, this corre- 
sponds to the transitions between the lowest and the 
highest adiabatic energy state: Si_>i = Sai^Au, an d 

Sn—*n 

— Sau^ai, where Ai are adiabatic states ordered 
by their energies (see Fig. 1). 

To compute these transition probabilities, I consider 
the n-component state vector ip(t) at large times t — + 
±00. At these times, the Hamiltonian is almost diago- 
nal, and the transitions between different adiabatic states 
are negligible. Therefore, the state i'th component ipi(t) 
evolves as 



1pi(t — > -OO) = CliC 

ipi(t — > +00) = bie 



-i<t>i{t) 

-iMt) 



(6) 



Let us consider the phase <j>i(i) = &rgipi(t). In the zero- 
order approximation, when one completely neglects the 
off-diagonal part A, the phase is entirely due to vc 



(*) 
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(7) 



To find the phase <j>i more accurately, one notes that in 
the adiabatic approximation this phase is actually given 
by an integral of a 'frozen eigenvalue' of ei(t) the Hamil- 
tonian (1): 



1 

i(t) = J £i(t')dt' . 



(8) 



At large t, vt ^> A, the eigenvalue Ei{t) can be found 
from perturbation theory in A. To the second order in 
A, it is given by 



e\ 2 \t) = v l t + A u +J2j ± 



{Vi - Vj )t 



(9) 



Therefore, at large t, the phase (fn{t) can expressed as 



M*) = ^f" + A u t + atlnt + 0{A 3 ) , (10) 



where 



(11) 



The key observation underlying our approach is that the 
terms retained in the expansion (10) are the only terms 
that grow indefinitely at t — > 00. The higher-order terms 
in A, written symbolically as 0(A 3 ), remain finite at 
large t. 

The asymptotic form of ipi(t), defined by the expres- 
sions (6), (10) and (11), allows one to compute the tran- 
sition amplitudes Su and S nn as follows [12]. Let us 
consider the solution (6) for arbitrary complex times t, 
allowing the phase variable (f>i(t) to be complex. Suppose 
that initially the system is in fc'th state, so that a% = 
for all i 7^ k, and a k — 1. At t — > +00, the system is, in 
general, described by all bi ^ 0. One can try to continue 
the solution at t — > +00 to t — > —00 along one of the two 
large half-circles, t = Re l9 : 

(a) in the upper half-plane, 9 £ [0, 7r], and 

(b) in the lower half-plane, 9 G [0, — tt]. 

The contour radius R is assumed to be very large. 
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FIG. 2. The contours in the complex time plane used for 
analytic continuation and asymptotics matching of the solu- 
tion for the i = 1, n components. 

When the state ip(t) is continued analytically along 
such a contour, its different components will behave dif- 
ferently: some of them increasing, while others decreas- 
ing as a function of \t\. Indeed, the leading term in <j>i{t) 
gives for the i'th component 



' ViR 2 , nn iuiR 2 
ipi oc exp ( — ^ — sin 29 H — cos 29 



(12) 



In matching the asymptotic expansions only the most 
rapidly increasing terms the wave function must be 



2 



retained. (For a discussion of this issue in depth, 
see [11,12].) For the contour (a) the leading term is given 
by the i = n component, while for (b) it is the i = 1 com- 
ponent. Because of that, only n'th and l'st components 
of the state ip(t — > — oo) can be determined in this way. 
Without loss of generality, I will consider the i = 1 case, 
i.e., perform continuation along the contour (b). 

The wave function after continuing along the contour 
can be found by changing the argument of t as 



t -> te 



(13) 



in Eq. (10). The first and the second term contribute 
only to the phase of ipi(t), while the third term changes 
its modulus: 



-oo) = b\e 



(14) 



On the other hand, we assume that this state evolved 
from the state which initially had only one nonzero com- 
ponent, |^>i(t)| = 1, fpj^i — 0. The time dependence of 
this component, ipi(t) oc e"'^W, with unit modulus and 
the phase given Eq. (10), yields 

tpi (t — > — oo) = exp (— l^it — iAnt — ia.\ hit) . (15) 

Comparing these two expressions, we obtain 

b\ = exp (wai) , (16) 

which gives the transition amplitude of the form 



iSi_>i = bi = cxp(7rai) = exp [ — 7r — 



• (17) 



Similarly, for the i = n component, after going through 
analytical continuation over the contour (a), we obtain 



S n ^n = CXp(-7ra„) = CXp -7T ^ 



I A, 



• (18) 



The resulting expressions for the transition amplitude, 
Eqs. (17), (18), bear similarity with the Landau-Zener for- 
mula and, in fact, coincide with it in the n = 2 case. Just 
like in the LZ problem, the dependence on the coupling 
matrix A and on the 'rapidity' v is such that S —> in 
the adiabatic limit of slow evolution (small v) , and 5 — > 1 
when the evolution is fast (large v). 

Eqs. (17) and (18) coincide with the transition ampli- 
tude found by Brundobler and Elser [10] in the sequential 



pairwise LZ level crossing approximation, and thus prove 
their conjecture. 

The existence of an analytic result for two transition 
amplitudes may indicate that the entire S'-matrix could 
be obtained analytically. It does not seem, however, that 
a trivial generalization of the above approach will be suf- 
ficient for that. We were not able to perform analytic 
continuation keeping track of the subleading terms, which 
probably indicates that this is not an optimal method of 
treating such a problem. 

It is of course not unconceivable, in principle, that only 
two transition amplitudes out of total n 2 can be obtained 
in a closed form, while others cannot. However, to the 
best of the author's knowledge, there are no other exam- 
ples known where only a subset of the 5-matrix is cal- 
culable, and thus such a possibility does not seem likely. 
The author hopes that this paper will stimulate further 
work towards a more complete understanding of multi- 
level Landau-Zener transitions. 
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